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NUMERICAL  SOLUTION  OF  FLOOD  rREDICTION  aI'^D 
KIVSR  ^:^GuLATIOF  PROBLEMS 

J,  J.    Stoker 

I.   Derivation  of  Basic  Theory  and  Formnlatlon 
of  Kuriieri c a  1  Methods  of  Attack 

1 ,   Introduction 

The  Institute  for  Mathematics  and  Mechanics,  llex-r   York 
University,  has  undertaken  the  task  of  investigatin'T  the 
possibility  of  using  modern  calculating  equipment  for  solving 
variolas  problems  conceminr:  wave  motions  In  rivers:   the 
problem  of  predicting  floods  from  the  observed  and  estiraatod 
run-off  from  precipitation,  the  problem  of  flood  stages  at  a 
jij.nction  of  two  large  ri.vers,  and  the  problem  of  predicting 
the  effects  due  to  operation  of  a  dam  arc  tyoical  examples. 
This  work  is  being  done  under  a  contract  with  the  Corps  of 
?;n.^'lneors  of  the  U.S.  Army  in  close  cooperation  with  the 
Ohio  River  Division  in  Cincinnati. 

The  Ohio  River  Division  already  has  a  satisfactory 
procedure  for  predicti.-iJ';;  floods  in  the  Ohio  on  the  basis  of 
the  observed,  or  estimated;  run-off  and  the  observed  flovj 
into  thu  Ohio  from  its  tributaries;  the  procoduj?o  is  called 
flood-routiiag.   The  flood-routing  procedure  can  bo  deduced 
as  an  approximation  to  the  solution  of  the  basic  differential 
equations  for  flow  in  opcui  char-ncls  (cf.  the  article  by 
B.  R.  Gilcrost  in  the  book  by  Rouse  [  6  ] '' " '^  )  j  but  it  makes  no 


^  ■  Square  brackets  rofcr  to  the  bibliography  at  the  end  of 
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direct  u/se  of  the-  cliff  or-cntial  oquat  ions .   Hoxrovor,  thcj  iloocl- 
roiitlng  proc;;chiri.;  In  qr'cr-i.ion  sGOrns  not  to  give  entirely 
sn.tisf actor:/  rosrilts  J.a  coses  other  than  tb.cit  of  ■:.lcterinining 
the  progress  of  a  flood  r'lcira   a  long  river--for  exanplc,  the 
eroblom  of  what  happens  at  a  junction,  such  as  that  of  the 
Ohio  and  Mississippi,  or  the  problem  of  calculatin'^;  the 
transient  effects  resulting';  from  regulation  of  a  dam,  such  as 
the  Kentucky  dam  at  the  mouth  of  the  Tennessee  river,  seem 
to  be  difficult  to  treat  by  muthods  that  are  modifications  of 
the  more  or  less  standard  flood-roiiting  proccdur-es.   Svcn  for 
the  Ohio,  the  usual  procedva^o  fails  occasionally  to  yield  the 
observed  river  stages,   '-in  thv.:  other  han^"! ,  the  basic 
differential  equations  fo]  flov/  in  open  channels  arc  ax-^plicablt 
in  all  cases  ancl  can  be  used  to  solve  the  problems  once  the 
appropriate  data  describing  the  river  and  the  appror'riate 
liiitial  and  boundary  coii'-hitions  are  knoi-ni.   The  present 
I'eoort  has  as  its  objocte   1)  the  derivation  of  -che 
differential  equations  and  a  discussion  of  a  nuj.iber  of 
properties  of  their  scl'.xtloris  which  arc  of  basic  importance 
for  a  proper  understanding  of  the  methods  to  b^.  used  to  solve 
the  concrete  probloma  TiUTiCrically ,  and  2)  a  discussion  of  the 
nuraorical  methods  to  bu  vscd  in  solving  the  problems. 

The  idea  of  using  the  differential  equations  di.i''octly 
as  a  means  of  treating  problems  of  flcv;  in  open  chaianols  is 
not  at  all  ncv;.   In  fact,  it  goes  back  to  Massau  [''.|.]  as  long 
ago  as  1889.   Since  then  the  idea  has  been  taken  up  by  many 
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othei's  (mostly  in  ip;nor^-ncc  of  the  woi'k  of  Mas3au)--ror 
example,  by  FroisvjGrk  [p] ,    ThoriiaK  [8],  Stoker  [7],   Thoinas, 
in  particular,  attached  the  flood-routing  problem  in  his 
noteworthy  and  pioneering  paper  [3]  and  oiitlinod  a  n.ui:.icrical 
procedur'j  for  its  3oli;.tio2i  based  on  the  idea  of  ur^inp"  the 
inethod  of  finite  differences.   However,  his  method  is  very 
laboi^ious  to  apply  and  X'/ould  .ilso  not  necessarily  fiu?nish  a 
f'ood  approximation  to  the  desired  solution  even  if  a  largo 
number  of  divirions  of  the  river  into  sections  v/crc  to  bo 
taken.   In  general,  the  araount  of  numerical  work  to  bo  done 
looked  too  formidable  for  practical  purposes  until  rather 
recently. 

During  and  since  the  late  war  new  devo  1  opine nt.«  have 
taken  place  which  make  th-e  idea  of  tackling  flood  prediction 
and  other  similar  proble::s  by  nijimcrical  solution  of  the 
relevant  differential  ccpAations  quite  tempting.   There  have 
been,  in  fact,  dovelopmoiits  in  two  different  directions,  both 
I'^otivatod  by  the  desire  to  solve  difficu.lt  problcras  in 
compressible  gas  dynamics:   1)  development  of  appropriate 
nuiaorical  proced-urcs  —  tor    the  most  part  iiiethods  using  finite 
diff erenc.js--for  solving  the  differential  eaua.tions,  and 
2)  dcvclopm.ont  of  comioTiting  machines  of  widely  varj'"ing 
characteristics  suitable  for  carrying  out  the  numerical 
calculations.   It  happens  that  the  differontial  equations  for 
flood  control  problems  arc  of  the  sam.e  type  as  those  for 
compressible  gas  dynamics,  and  consequently  the  cxporionco 
and  calculating  equipment  developed  for  gas  dynam.ics  ;^;roblems 
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can  be  used,  or  Guita':.0_y  r.io'lif  icd ,  for  solving  flood  control 
problems.   The  basic  gonoral  objocbivo  of  tho  work  of  the 
IMM  in  this  field  in  to  analyze  those  possibilities  and  tost 
thon  by  coding  concrete  problems  for  sviitublc  calculating 
machines  and  comparing  the  results  ^^7ith  observations. 

In  section  2  a  brief  derivation  of  the  well-known 
differential  equations  governing  flovr  in  open  chann;:lc  will 
be  given.   In  section  3  various  interesting  and  basic  general 
proriertios  of  tho  solutions  of  the  differential  equations 
are  discussed  on  the  basis  of  the  theory  of  characteristics. 
It  is  possible  to  make  use  of  the  characteristic  form  of  the 
d5.ffercntial  equations  to  solve  them  nuuiiCirically,  out  \-rc 
sh.all  not  proceed  in  siach  a  fashion.   Never  tho  less  an  under- 
standing of  the  theory  ox  the  equations  based  on  the  use  of 
the  characteristic  form  is  vital  for  th';  proper  vu'.tder standing 
and  use  of  tho  numerical  r-et]T.ods  we  propose  to  use.   In 
section  I4.  the  nurnorical  '-lethods  we  employ  are  doscribod; 
thoy  are  methods  in  wi-icli  Wio  fi^'le,  of  integration  in  tho 
x,t-planc  is  covered  bj?-  a  rectangular  n^^t  of  points,  the 
derivatives  in  the  differential  equations  ar-^  replaced  by 
difference  quotients,  luid  the  approximate  solutions  are 
foTond  by  solving  ].incar  oqviations  for  the  values  of  t;io 
desired  quantities  ( ossentiailj^  the  depth  of  tho  water  and 
its  flow  speed)  at  tho  net  points. 


5. 

Z ,   Brief  derivation  or  tlic  f'if  f  orcntlal  oquationa 

The  differential  cqr'.ations  T^?hich  form  the  basis  for 
our  subseqi.vent  dincussion  r.ro  uoll  known  and  many  derivations 
of  fchom  have  been  given  (of.,  for  example,  the  article  by 
Gilcrcst  [6]).   For  the  sake  of  completeness,  ajid  to  fix 
notations,  we  give  a  brief  derivation. 

The  theory  is  one-diimcnsional ,  i.e.  the  actual  flow 
in  the  channel  is  assumed  to  be  well  approxim.atod  by  a  flow 
T:ith  uniform  velocity  over  each  cross-section,  and  the  free 
surface  is  taken  to  bo  a  level  line  in  each  cross-section. 
The  ri-ver  is  assumed  also  to  be  straight  enough  that  its 
course  can  bo  developed  into  a  straight  line  without  causing 
serious  errors  in  the  flow.   The  flow  velocity  is  denoted 
by  V,  the  depth  of  the  stream  by  y,  and  these  quantities  are 
functions  of  the  distance  :c  down  the  stream  and  of  the 
time  t.  (cf.  Fig.  2.1).   Vh.e  vertical  coordinates  of  the 
bottom  and  of  the  free  siix-facc  of  the  ftroam,  as  measured 
from  the  horizontal  a:cis  x,  are  denoted  by  h(x,t)  and  z(x). 
The  breadth  of  the  free  surface  at  any  section  of  the  stream 
is  denoted  by  B. 
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Basic   wotatlons 


7. 

Tho  ciffcrontial  oqva.tions  govcrninp;  the  flew  arc 
expressions  of  the  laws  of  conservation  of  mass  and  rj:;raontum. 
In  deriving  thorn  the  follov/ing  assumptions,  in  addition  to 
those  mentioned  above,  are  iiiade^'':   1)  the  pressiiro  in  tho 
water  obeys  the  hydrostatic  pressure  law,  2)  the  slope  of  the 
bed  of  the  river  is  small,  3)  the  effects  of  friction  and 
turbulence  can  bo  accountec"'  for  tiirough  the  Introduction  of  a 
resistance  force  depending  on  the  sauare  of  the  velocity  v 
and  also,  in  a  certain  T/ay  to  be  specified,  on  tho  depth  y. 

Wo  first  derive  tl-c  so-called  equation  of  continuity 
froin  tho  fact  that  the  riass  oAAx  included  in  a  layer  of 
x-70.ter  of  density  p,  thiclrncss  Ax,  and  cross- sect  ion  area  A, 
changes  in  its  flow  along  the  stream  only  througli  a  possible 
inflow  along  tho  banks  of  the  stream,  say  at  the  rate  pq  per 
unit  length  along  the  river:   here  wo  use  the  fact  that  tho 
particles  in  a  vertical  plane  section  remain  in  it  by  virtue 
of  on^  of  the  assiamptions  'iiado  above.   This  moans  that  the 
so-called  particle  derivative  -^(pAAx)has  the  value  pqAx, 
and  honco  that  the  equ.;.tion 

(2.1)  v„AxA  +  A;:^-i  +  vA  )  =qAx   holds   , 


These  assumptions  are  not  the  minimujn  number  necessary: 
for  example,  assijimptlon  1)  lias  as  a  consequence  the 
independence  of  the  velocity  on  the  vertical  coordinate  if 
that  wore  true  at  any  one  instant  (cf.  [?],  p.  lii.  and  p.  6l) 
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since  4i-  =   C,    —-x~  -   v.,lS'Z  roprosents    the    change    in   thickness 

r'  i'  -;;- ) 

of    the    slice,    and    th,    o^ciTLtion  —  is   .''efincd  bj^ 

fit 

(2.2)  ._L..  =  -4.  +  V -^ 

^        '  Ct       e>c         $x        ' 

with  V   the    velocity    oi    tlio    strcj.m.     Equation    (2.1)    c^.n  be    put 
In  the    forra 

(2.3;  (Av),.   +   A^    =   q 

Here  q  represents  th^  volume   of  the  flow  into  the  river 
over  its  banks  ■ecr  unit  length  and  tine.   (This  eouation 
coiild  also  be  dorived  ""ircctly  and  easily  by  calculating  the 
flov;  into  and  out  of  the  slice  of  thickness  Ax.)  i^'or  the 
special  case  of  a  rectarj.,-?ular  cliannel  of  constant  br'^adth  B, 
and  hence  such  that  A  =  By,    the  equation  of  continuity  is  at 
once  soon  to  be  expressible  as  follows: 

(2.b-)  ^x^'  "^  "^^x   "^  '^t  "  ^^^ 

The  equation  of  v.iotion  is  next  cerivcd  for  the  same  mass 
incliided  between  tvjo  vertical  planes  a  distance  Ax  apart. 
Tnc   momentum  of  such  a  .nass  is  given  by  oaAxv,  and  Newton's 
second  law  states  that 


-pF) 

^  '  If  ^(Xjt)  is  an;,'-  function  associated  with  a  particle,  i.e. 

attached  to  a  particle,  v-jhoso  motion  is  given  by  x=x(t),  it 
follows  that  §=Q^  g-.Q.  =  (^^v^)Q^  sine,  f  =  v,  the 

velocity  of  the  Tjarticle .   Thus,  since  Z\-X  =  x ,  -x  ,  with  x,  the 

+        —  + 

coordinates    of   the    vortical   planes    c'efining    our    slice,    it 

'^  A  X         "a 

follows  that  —^TT~  =  'J',-'   V   -^  v  Ax  since  x  is  one  of  the 

indopcndGnt  variables.   It  should  be  observed  that  the-  sub- 
scripts X  and  t  are  used  to  denote  partial  derj.vatives , 
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(2.5) 


4(o..Axv) 
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vjith  -^^  onso  more  th.'-;  r'^.rfciclo  doi-ivativo,  and  F  the   resultant 
of  all  forcer,  on  the,  nlicc  in  the  horfizontal  f.lroction.   The 
forcoE  acting  on  tho  clci.iont  in  question  arc  IndicatG;''.  in 
Pig,  2,2.   The  force  P  is  t}ie  resultant  pressure  nn   c.   vertical 
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cross-SGCtion,  tho  forcoG  r,  represent  the  pressure  in  the 
horizontal  plane  frorn  t-ic  river  banlis,  G  is  tho  xicight  of  the 
element  J  F^  tho  friction  force,  and  N  the  prossiirc  normal  to 
tho  river  bottom.   3ecau3C  of  tho  hydrostatic  prcssiij:"© 
assumption  we  have  (cf,  I'^'ig.  2.1) 

(2.6)  ^  ""  f  PS'>^-  1''°^^^'*^^"  ^  ' 

and  hence  the  quantity  P.,  is  given  by 

(2.7)  P^  =  r  pgy^b(;c,  r|)dr^  +  f  ps(  y  -  rj  )b_/,x,  t<|_)  d  »! 

^0  ^0 

J 

The  second  term  on  tho  i-'irh,t  can  be  interpreted  as  the  x- 
component  of  tho  pressure  P,, ,  and  in  thus  balanced  by  this 
pressure.   Since  -^(pAOx)  =  paLSx   is  tho  continuity 
condition,  it  follows  that  we  have  from  (2.5)  the  equation 

(2.8)  pAZAx(v,  +vv  )  +  pqAx'V  =  Z\  xN  sin  (j)  -  Ax?^  coscb-P  Ax. 

1^      X  i  X 

Since  the  slope  angle  (j)  Is  small,  wo  can  replace  siii  <j?  by 
"T^  =  S,  the  slope  of  bhe  river  bed,  and  cos  (j)  by  1.   In 
addition,  N  can  be  replacod  b3'  G--again  because  ij)  is  small 
--and  hence 


(2.9)  Axil  =  AxApi 


For   the   resistance   force  F^  we   use   Manning's   formula: 


IV    ■ ,  (  , 
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(2.10)  F.  =  ^^^^^        , 

which  states  that  the  posifitance  is  opposite  in  direction  to 
the  velocity,  but  is  proportional  to  its  square.   The  constant 
Y  is  a  roughness  coefficient,  and  R  is  the  hydraulic  radius, 
i.e.,  the  cross-section  area  divided  by  the  wetted  perimeter 
(cf.  [6]).   Thus  we  obtain,  finally,  the  equation  of  motion: 

(2.11)  V,  +  vv_  H-S'v  =  Sg  -  gy  - -^  F^:, 

=  Sg  -  S^g  -  gy^        , 
X\rith   S^   the    so-called   friction   slope    defined  hj 

(2.12)  S..   =  ~  F„ 

X  pgA      f 

The  differential  equations  (2.3)  and  (2.11)  arc  the 
differential  equations  which  are  basic  for  the  v/orl:  to  follow. 
Of  course,  the  cross-section  area  A  must  be  regarded,  as  a 
known  function  of  y  that  is  fixed  by  the  topography  of  the 
river  valley;  thus  (2.3)  coiild  also  be  written  in  the  form 

(  2  .  13  )  ^V^r^  "'"  '''-'^  ^  ^^v  "^  ^.r^-  =  ^ 

J  -/».      -(i.        j\.  y      v 

and  (2,11)  and  (2.13)  then  arc  a  pair  of  equations  to 
determine  the  two  quantities  y(x,t)  and  v(x,t). 

For  any  given  river  it  is  thus  necessary  to  have  data 
available  capable  of  yieldLing  the  cross-section  area  A,  the 
hydraulic  radius  R,  and  the  roughness  coefficient  y  of  the  river 
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as  functions  of  x  and  j,    aiid  of  tho  slope  S  of  its  bod  as  a 
fujiction  of  X  in  order  to  have  the  coefficients  in  the 
differential  equations  (2.j,l)  and  (2.13)  defined.   TVirce  of 
those  quantities  are  purely  geometrical  in  character  and 
could  in  principle  bo  dotcriii.ined  from  an  acc^xratc  contour 
y:\e.p   of  the  river  valley,  but  tho  determination  of  the 
rou,ip;iiness  coefficient  y  '^^   course  requires  mcasiiroiiionts  of 
actual  flov;s  in  the  river  for  its  determination.   (These 
remarks  imply  that  the  river  does  not  change  its  course  fr 
time  to  timc--as  happens,  for  example j  in  the  Misso-'jj?i--30 
that  the  coefficients  of  the  differential  equations  arc 
independent  of  the  time . ) 


3.   Tho  method  of  characteristics 

It  has  already  been  stated  that  many  important 
properties  of  tho  solutions  of  our  differential  equations  and 
ma.ny  valuable  insights  with  regard  to  the  integration  th^^ory 
associated  with  thorn  can  be  gained  through  the  use  of  v/nat  is 
called  the  method  of  characteristics.   In  this  section  we  give 
an  explanation  of  this  method,  but  since  vro  do  not  make  direct 
use  of  it  later  on  in  solving  problems  nurrierioally  v;c  shall 
base  cur  explanation  on  a  s-.ecial  case  of  the  differential 
equations,  although  the  results  and  observations  are  applicable 
in  the  most  general  case.   The  special  case  in  question  is 
that  of  a  river  of  constant  rectangular  section  and  uniform 
slope,  with  no  flow  inLo  the  river  from  the  banks  (i.'  ,  q  =  0 
in  (?.ll)  and  (2.13)).   In  this  caso  the  dlff crontic:i 
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equations  (2.11)  c.nd  (2.13),  with  A  =  5y,  can   be  written  as 
follows : 

(3.1)  V  y  +vy.^.  +y,  =  0 

(3.2)  v^  +  vv,.  +  gy,^  +  E  =  0 

!\V.  have  introduced  the  sy-mbol  E  for  the  external  forces  per 
u:'.\it  mass: 

(3.3)  E  =  -gS  +gS^   ,   S  =  const. 

The  term  E  differs  frora  the  others  in  that  it  contains  no 
derivatives  of  y  or  v. 

The  theory  of  chai^acteri sties  for  those  equations  can 
be  approached  very  directly^    in  the  present  special  case  by 
introducing  a  new  quantity  c  to  replace  y,  as  follown: 

This  qiaantity  has  groat  physical  significance,  einco  it 
ropresents--as  we  shall  see--the  propagation  3]~/0cd  of  small 
disturbances  in  the  river.   From  (3.-!-)  'w^  obtain  at  once  the 
relations 


(3.5)  2cc   =  ^y,^   ,    2cc   =  gy^ 


and  the  differential  equations  (3.1)  '^-nd  (3.2)  take  the  form 


T-^n ■  ■ 

*-  '  For  a  trcatmunt  v/nich  shows  quite  generally  how  to  arrive 
at  the  formulation  of  the  characteristic  equations,  sec  [1], 
Ch.  2. 


1^. 


"  2  c  c  +  V  ,  +  V  V  +  E  =  0   , 


(3.6) 


\cv  +2vc  + 


2c.  =  0 


Tliese  equations  are  next  added,  then  subtracted,  to  obtain 
the  following  equivalent  ^:'a.ir  of  equations: 


. ,.  •..■-■,  V  +  E  -    0 


(3.7) 


(--2{(-c+v)^.-^}cHj-c-^v)^  +  ^}v+E 


Ue  observe  that  the  derivatives  in  these  equations  now  have 
the  form  of  directional  derivatlves--indced,  to  achieve  that 
vras  the  purpose  of  the  trajisformation--so  that  c  aiid  v  in  the 
first  equation,  for  exara'olc ,  are  both  subject  to  f'xo    operator 
[O^v)—  +  —p,    which  means  (cf,  the  footnote  in  the  ■'•jreced  tng 
section)  that  these  func'ricns  are  differentiated  along  curves 
in  the  x.,t-Y>lane  which  satisfy  the  differential  equati(;n 

dX 

^,-  =  c+v.   In  similar  fashion,  the  functions  c  and  v  in  the 
second  equation  are  bot.j.  subject  to  differentiation  along 
cu.rvo3  satisfying  the  difforontial  eouation  -r^r  =  -c+v. 

It  is  entirely  feac'ible  to  develop  the  integration  theory 
of  equations  (3«7)  quite  .~ene rally  on  tjie  basis  of  these 
observations  (as  is  done,  Tor  exaraple,  in  [1],  Cii.  2),  but  it 
is  sirapler,  and  leads  to  the  same  general  results,  to  describe 
it  for  the  special  case  in  v;hich  the  resistance  force  F .,   is 
neglected  so  that  the  quantity  S  in  (3.7)  is  a  constan.t  (see 
(3o/).   In  this  case  the  equations  (3.7)  can  be  vn-'itton  in 
the  form 


■  'F 
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(  ^  + (v+c)^-.-!  (\+2c  +Et)  =  0 
(3.7)i 

as  one    can  readily  verify.   But  the  interpretation  of  the 
operations  defined  in  (3»7)-]  'tias    just  been  discussed:   The 
relations  state  that  the  f">inctions  (v±2c+Et)  are  constant 
for  a  point  moving  through  the  fluid  with  the  velocity  (v±c), 
or,  as  V7e  raay  also  put  it,  for  a  point  vfnose  motion  in  the 
:c,t-plane  is  characterizod  by  the  ordinary  dif forcntial 
equations  -yr-   =  v±c.   That  is,  wo  have  the  follox;ing  situation 
in  the  x,t-plane:   Thoro  arc  two  sets  of  curves,  C-,  and  Cp 
called  characteristics,  u-iich  are  the  solution  curves  of  the 
ordinary  differential  equations 


^1-    i:^  "  ^'"-^    '    -''' 


r 


-         h.-^:.  :r  v-C 


.nd  V70  have  the  relations 

v+2c  +Et  =  1:^  =  const,  alon^  a  curve  C-i  an' 
(3.9) 

v-2c  +  iJt  =  1:.-,  =  const,  alon?-;  a  curve  Co 

Of  course  the  constants  k^  and  kp  will  be  different  or. 
different  curves  in  general.   It  should  also  be  observed  that 
the  tv;o  families  of  characteristics  determined  by  (3.3)  are 
really  distinct  because  of  the  fact  that  c  =  /gy  7^  u  since  wc 
suppose  that  y  >  0,  i.e.  tnab  the  water  surface  never  touches 
the  bottom. 


''■  ,:■■ 


By  reversing  the  abovo  procedure  it  can  bo  coon  rather 
easily  that  the  systei-i  ci  relations  (3.8)  and  (3,9)  is 
coinpietely  ecuivalont  to  bho  system  of  equations  {3*^''i    fo^  "the 
case  of  constant  bottom  sIolo  and  zero  resistance,  so  that  a 
solution  of  either  syr.tcm  yields  a  solution  of  the  othon^ ,   In 
fact,  if  we  set  f(x,t)  -   v-!-2c  +  Et  and  observe  that 
f(x,t)  =  k^  =  const,  along  any  curve  x  =  x(t)  for  vrhich 
^  =  v+c  it  follows  that  along  such  curves 


(3.10)  f^^f^H.  f^^(v+c)f^^  =  0    . 

In  the    same  way  the   fi-.nction  g(x,t)    =   v-2c  +  St    satisfies   the 
relation 

(3.11)  g^  ^•  '',v~c)g^   =   0 

dx 
along  the  curves  for  ijhicli  --rT-  =  v-c.   Thus  vjherevcr  the  cm''ve 

families  C-,  and  Cp  cover  ':i-'.o  x,t-planc  the  relations  (3.10) 
and  (3.11 1  hold.   If  nov;  equations  (3,10)  and  (3.11)  ^-^o   added 
and  the  definitions  of  f(::^t)  and  g(x,t)  are  recalled  it  is 
readily  seen  thc^t  the  firr:t  oT   equations  (3»t>)  results.   By 
subtracting  (3.11)  froiri  (3,10)  thu  second  of  aquations  (3.6) 
is  obtained.   In  other  worlds,  any  functions  v  and  c  '''hich 
satisfy  the  relations  (3.'-)  ^-!^<^  (3.9)  will  also  satisfy  (3,6) 
and  the  tx-ro  systems  of  equations  are  therefore  now  seen  to  be 
c crap lo  t  c  ly  e  qu i va  1  ont . 

As  wo  Xifould  expect  on  physical  grounds,  a  solution  of 
the  original  dynamical  equations  (3.6)  could  bo  shoxrn  to  be 
uniquely  determined  vih.cn   appropriate  initial  conditions 


rr'^0 


17. 

(for  t  =  0,  say)  and  coT,indJ.ry  conditions  arc  precoribed;  it 
fol3_ov7s  that  my  soltvfclons  c±  (3,8)  and  (3.9)  arc  also  unique 
when  such  conditions  arc  v.rc-scribod  since  i\fo  kno;;  that  the  tvjo 
systems  of  equations  arc  oquivalent. 

At  first  sight  one  night  he  inclined  to  rcgai-^d  the 
rc?.a.tions  (3.8)  and  (3  •9)  as  more  complicated  than  the  original 
pair  of  partial  differential  eqiiations,  particularly  since 
the  right  hand  sides  of  (3«8)  S-^w  not  known  and  hence  the 
characteristic  curves  are  also  not  known.   Nevertheless,  the 
formulation  in  terms  of  the  characteristics  is  quite  i.^seful 
in  studying  properties  o.C  tiie  solutions  and  also  in  studj^ing 
questions  referring  to  the  appropriateness  of  various  boundary 
and  initial  conditions.   It  is  useful  to  begin  by  describing 
briefly  a  method  of  dotorjiiining  the  characteristics  by 
successive  approxiriiations  vjhich  at  the  same  time  mak'os  possible 
a  niraber  of  useful  observations  and  interpretations  regarding 
the  role  of  the  characteristics  in  general.   Let  us  consider 
a  problem  in  which  the  values  of  the  velocity  v  and  t"io  surface 
elevation  y  (or,  v;hat  a:;io-u.nts  to  the  same  uhmg,  the  ijave 
speed  c  =  fi^)arc   proscribed  for  all  values  of  x  when  t  =  0. 
k'e  lAiish  to  calculate  the  solution  for  t  >  0  by  determining  v 
and  c  through  use  of  (3.n)  and  (3.9)  and  the  given  initial 
conditions.   We  have  for  t  =  0  the  conditions 


(3.12) 
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in  which  v{x)    and  'c{x)    arc  given  functions.   V/c  can  approximate 
the  values  of  v  and  c  for  small  values  of  t  as  follows: 
consider  a  series  of  points  on   the  x-axis  (cf.  Pip;.  3»1)  a  small 
distance  Ax  apart.   At  all  of  these  points  the  values  of  v 
and  c  are  known  from  (3.12).   Ccnsoquontl  j'  the  slopes  of  the 
characteristics  G-,  and  C-^  at  tVieso  points  are  known  from  (3.8). 
From  the  points  1,  2,  3,  U  straight  line  segments  with  tb.cse 
slopes  are  drawn  until   they  intersect  at  points  5.-  6»  and  ?, 
and  if A^is  chosen  sufficiently  small  it  is  reasonable  to  expect 
that  the  positions  of  these  points  will  be  good  approximations 
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Fig.  3.1 

to  the  intersections  of  the  characteristics  issuing  from  the 
points  1,  2,  3,  1^.   The  values  of  both  x  and  t  at  points  5,  ^, 
and  7  are  now  known- -they  can  be  determined  graphically  for 
example --and  through  the  use  of  (3.9)  and  the  initial 
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conditions  we  can  also  determine  the  approximate  vo.liies  of  v 
and  c  at  these  points.   For  this  purpose  v:ti  observe  that  along 
any  -oarticular  segment  issuing  rrorti  the  points  1,  2,  3  or  \\ 
the  values  of  v-i-2c  +  ZIt  and  v-2c  +  Et  are  known  constants  since 
the  values  of  v  and  c  are  fixed  by  (3.12)  for  t  -   0;  henoc  we 
have 

r  Along  ^'■-,  :  v+2c  +5t  =  v  +  2o   ,   and 

I  Along  C^:  v-2c+St  =  v  -  2 c   . 

At  the  points  5?  6,  and  7  '.'c  icnow  the  values  of  t  ane  hence 
(3.13)  furnishes  two  independent  linear  equations  f.ov    the 
dctcrraination  of  the  values  of  v  and  c  at  each  of  those  points 
Cnco  V  and  c  are  l-cnown  at  points  $,    6,  and  7  the  slo-jcs  of  the 
characteristics  issuing  froiu  these  points  can  be  deterir.ined 
from  (3.8)  and  the  entire  j^rocess  can  be  carried  ov'c    once  norc 
to  yield  rhe  additional  vjoints  S  and  9  and  the  approximate 
values  of  v  and  c  at  theso  points.  Ln   this  way  vre  oan 
approximate  the  values  of  v  ana  c  at  the  points  of  a  net  over 
a  certain  region  of  the  Xjt-plane,  and  can  then  obtain 
approximate  values  for  v  r::c\sj_   c  at  any  points  in  the  sarac 
region  either  by  interpolation  or  by  rofirxing  the  net  insido 
the  region.   It  is  quite  plausible  and  could  be  proved 
raathcmatlcally  that  the  above  process  would  converge  as 
A>^  — ^  0  to  bho  unique  solntion  of  (3«S)  and  (3  •9),  -nd  hence 
also  of  (3.6),  ccrrosponding  to  the  given  initial  conditions 
for  sufficiently  small  values  of  t  (i.e.  for  a  region  of  the 
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x,t-plane  not  too  far  from  the  x-axis)  provided  that  the 
proscribed  initial  val^ios  of  v  and  c  are  sufficiently  regular 
fimctions  of  x. 

It  should  bo    clear  that  once  the  characteristics  are 
laio\-jn   the  values  of  v  and  c  for  all  points  of  the  x,t-plane 
covered  by  then  are  also  Icioi'm,  since  the  constants  k^  and  kp 
in  (3 •9)  are  known  on  each  characteristic  through  the  initial 
data  and  hence  the  values  of  v  and  c  for  any  point  (x,t)  can 
be  calcvilatcd  by  solving  the  linear  equations  (3 •9)  for  the 
characteristics  through  that  point.   This  stateraent  of  course 
implies  that  each  one  of  the  two  families  of  characteristics 
covers  a  region  of  the  :c,t-plane  simply  and  that  no  tivo 
members  of  different  fariilies  are  tangent  to  each  other  —  in 
other  vjords  it  is  implied  tliat  the  two  families  of  character- 
istics form  a  non- singular  curvilinear  coordinate  system  over 
the  region  of  the  x,t-plane  in  question. 

The  method  of  finite  differences  used  above  to  determine 

the  characteristics  can  be  interpreted  in  such  a  way  as  to 

throiv  light  on  the  physical  properties  of  the  solution. 

Consider  the  point  10  of  Fig.  3«1  for  example.   We  recall  that 

the  approximate  values  v^  ,  and  c,  ^  of  v  and  c  at  rjoint  10  were 

i  o       i  w 

obtained  through  makin'^:  use  of  the  initial  values  of  v  and  c 
at  points  1,  2,  3,  [j.  on  the  x-axis  only,  and  furthermore  that 
the  values  v,„  and  c-,  ^  required  the  use  of  points  confined 
solely  to  the  region  v;ithin  the  approximate  characteristics 
Joining  point  10  with  points  1  and  '+.   Since  the  finite 
difference  scheme  outlined  above  converges  as  Ax  — *«  0  to 
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yield  the  exact  characteristics  it  is  clear  that  i-7C  can  make 
the  follovjing  important  statement:   the  values  of  v  and  c  at 
any  point  F(x,t)  within  the  region  of  oxistonco  of  the  solution 
arc  determined  solely  by  the  initial  values  prescribed  on  the 
segment  of  the  x-axis  x-j-hich  is  subtended  by  the  two 
char'actcristics  issuing  from  P.   In  addition,  the  two 
characteristics  issuing  from  P  are  also  determined  solely  by 
the  initial  values  on  the  segment  subtondcd  by  them.   Such  a 
segment  of  the  x-axis  is  often  called  the  domain  of  dependence 
of  the  point  p.   Correspondingly  wc  may  define  the  range  of 
influence  of  a  point  Q  on   the  x-axis  as  the  region  of  the 
x,t-plane  in  which  the  values  of  v  and  c  are  inflticnccd  by  the 
initial  values  assigned  '^o  y^oint  Q.   In  Fig.  3»^  ^'f'-^  indicate 
those  two  regions. 
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VJe  arc  now  in  c.  position  to  understand  the  role  of  the 
characteristics  as  ciirvcs  along  which  discontinuities  in  the 
first  and  higher  doriyatlvos  of  the  initial  data  are  propagated, 
since  it  is  reasonable  to  expect  (and  could  be  proved)  that 
those  points  P  vjhosc  domains  of  dependence  do  not  contain  such 
discontinuities  are  points  at  which  the  solutions  v  and  c  also 
have  continuous  derivatives.   On  the  other  hand,  it  could,  bo 
shoTim  that  a  discontinuity  in  the  initial  data  at  a  certain 
point  does  not  in  general  die  out  along  the  characteristic 
issuing  from  that  point.   Such  a  discontinuity  (or  distiorbancc 
in  the  water)  therefore  spreads  in  both  directions  over  the 
s-.ii''face  of  the  water  vjith  the  speed  v+c  in  one  direction  and 
v-c  in  the  other  in  view  of  the  interpretation  given  to  the 
chara.cteristics  through  the  relations  (3  •?)-]•   Sixice  v  is  the 
velocity  of  the  water  particles  we  see  that  c  represents  quite 
generally  the  speed  at  wliicn  a  discontinuitj'"  in  a  derivative 
of  the  initial  data  propagates  relative  to  the  moving  x^rater  . 
V7o  are  therefore  justified  in  referring  to  the  quantity  c  =  .J  gy 
as  the  ^^?ave  speed  or  propa.gation  speed. 

¥c  considered  above  a  problem  in  which  only  Initial 
conditions,  and  no  boundary  conditions,  were  prescribed.   In 
the  problems  wo  consider  later,  hovjevcr,  such  boundary 
conditions  vrill  occur  in  the  form  of  conditions  prescribed  at 
a  certain  fixed  point  of  the  river  in  terras  of  the  time: 
for  example,  the  height,  or  stage,  of  the  river  might  bo  given 
at  a  certain  station  as  a  function  of  the  time.   In  other  words, 
conditions  would  be  prescribed  not  only  along  the  x-a:cis  of 
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oijir  x,t-pl-anc,  but  also  o.long  the  t-axis  (in  gcncpal  only 
for  t  >  0)  for  a  certain  fixed  valuo  of  x.   The  method  of 
finite  differences  used  above  to  discuss  the  Initial  value 
problem  can  be  modified  in  an  obvious  x-ray  to  deal  i/ith  cases 
in  which  boundary  conditions  are  also  imposed.   In  doing  so, 
it  would  also  become  clear  just  what  kind  of  bomidary 
conditions  could  and  should  be  imposed.   For  example,  in  the 
great  majority  of  rivers--in  fact,  for  all  in  which  the  flow 
is  sub-critical;  i.e.  such  that  v  is  everywhere  less  than  the 
\-jo.vc    speed  Jgy  --it  is  possible  to  prescribe  only  onc_ 
condition  along  the  t-axis,  which  might  be  cither  the  velocity 
V  or  the  depth  y,  in  contrast  with  th.:^  necessity  to  impose 
two  con'^.itions  along  the  x-axls.   This  fact  would  become 
obvious  en  setting  up  the  finite  difference  scheme,  and 
examples  of  it  will  be  seen  in  later  reports. 

Finally,  it  should  be  stated  that  the  role  of  the 
characteristics,  and  also  the  m.cthod  of   finite  differences 
applied  to  them  could  bo  used  with  reforcnce  to  the  general 
case  of  the  characteristic  equaticriuS  as  embodied  in  tlac 
equations  (.3.7)  and  (3.e)  in  essentially  th^j    same  v;ay  as  was 
sketched  out  above  for  the  system  comprised  of  (3»5)  and 
(3 •9)  which  referred  to  a  special  case.   In  particular,  the 
role  of  the  characteristics  as  curves  along  x-ihich  small 
disturbances  propagate,  and  their  role  in  determining  the 
domain  of  dependence,  range  of  influence,  etc.  remain  the 
same . 
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I4..   Approximate  methods  to  be  used  in  c^ilculatin,^  solutions 
of  tho  dlfforcntli\?L  oquatlons 

It  has  alreadj'-  been  stated  that  whllo  the  f omiilatlon 
of  our  problems  by  the  riiothod  of  characteristics  is  most 
valuable  for  studying  many  questions  concerned  with  the 
solutions  of  the  differential  equations,  it  is  in  most  cases 
not  tho  best  formulation  to  use  for  the  purpose  of  calculating 
the  solutions  niimcrically .   Tnat  is   not  to  say  that  the 
device  of  replacing  derivatives  by  difference  quotients  is 
given  up,  but  rather  that  this  device  is  used  in  a  different 
manner.   The  basic  idea  is  to  operate  with  finite  differences 
by  using  a  fixed  rectangular  net  in  the  Xjt-planc,  in  contrast 
with  the  method  outlined  above,  in  which  tho  net  of  points 
in  the  x,t-planc  at  wliich  the  solution  is  to  be  approximated 
is  determined  only  gradrially  in  tho  course  of  the  computation. 
In  the  latter  procedure  it  is  thus  necessary  to  calciilate  not 
onlj   the  values  of  the  vjiknown  f\inctions  v  and  c,  but  also  the 
valines  of  the  coordinates  "yt   of  the  net  points  thomsolves, 
vjhcrcas  a  procedure  making  use  of  a  fixed  net  vjoule  require 
the  calculation  of  v  a.nd  c  only,  and  it  i>;Guld  also  have  the 
advantage  of  furnishing  these  values  at  a  convenient  sot  of 
points . 

However,  the  question  of  the  convergence  of  tho 
approximate  solution  to  the  exact  solution  v;hon  tho  mesh  width 
of  a  rectangular  net  is  made  to  approach  zero  is  more  delicate 
than  it  is  when  the  method  of  characteristics  is  used.   For 
example,  it  is  not  correct,  in  general,  to  choose  a  net  in 
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wliich  the  ratio  of  the  niosh  width  At  along  the  t-axic  and 
the  mesh  width  Z\x  along  the  x-axis  is  kept  constant:   such 
a  procedure  would  not  in  general  yield  approxiraations 
converging  to  the  solution  of  the  differential  equation 
"oroblcm.   The  reason  for  this  can  be  understood  with  reference 
to  one  of  the  basic  facts  about  the  solutions  of  the 
differential  equations  vrhich  was  brought  out  in  the  discussion 
of  the  preceding  section.   The  basic  fact  in  question  is  the 
existence  of  what  was  called  there  the  domain  of  dependence 
of  the  solution.   For  example ,  suppose  the  solution  x-rcre  to 
bo  approximated  at  the  points  of  the  net  of  Fig.  Li-.la  by- 
advancing  from  one  rov;  parallel  to  the  x-axis  to  the  next 
vo\-j   a  distance  At  from  It.   In  addition,  suppose  this  were 
to  be  done  hj   determining  the  approximate  values  of  v  and  c 
at  any  point  such  as  P  I'cf.  Fig.  Ij-.lb)  by  using  the  values 
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of  these  quantities  at  the  nearest  three  points  J,  1,  2  in 
the  next  line  below,  replacing  derivatives  in  tiic  tvjo 
differential  equations  by  difference  quotients,  end  then 
solving  the  resulting  algebraic  equations  for  the  tvro  ujnknowns 
V-  and  c^ .   It  soems  reasonable  to  su.ppose  that  such  a  scheme 
would  be  appropriate  only  if  P  were  in  the  triangular  region 
forricd  by  the  characteristics  drawn  from  points  0  an;'  2  to 
form  the  region  x-zithin  w'lich  the  solution  is  determined 
solely  by  the  data  given  on  the  segment  C-2:   otherv/ise  it 
seems  clear  that  the  values  at  additional  points  on  tiio 
X-axis  ought  to  be  utilized  since  our  basic  theory  tolls  us 
that  the  initial  data  at  some  of  them  vjould  indeed  influence 
the  solution.   On  the  etlier  hand,  the  characteristic  curves 
depend  themselves  upon  the  values  of  the  unlcnovm  fujicvions 
v  and  c  — their  slopes,  in  fact,  are  given  (cf.  (3»3))  by 
—■  ~  V  i  c  and  thus  the  ini^crval  l\t   must  be  choson  srar.ll 
enough  in  relation  to  a  fixed  choice  of  the  interval  l\x  so 
that  the  points  such  as  ir  •.:ill  fall  within  the  apprc'^riate 
domains  of  determinacy  relative  to  the  points  used  in 
calculating  the  solution  at  P.   In  other  words,  the  theory 
of  the  characteristics,  even  if  it  is  not  used  dir.;ctl3'-, 
comes   into  play  in  deciding  the  relative  values  of  /-LXt  and 
l\x   for  convergence.   (^'liis,  by  the  way,  is  the  reason  for 
the  statement  made  earlier  that  the  finite  difference  cchomo 
proposed  by  Thomas  [3]  mi/'ht  fail  to  converge  no  matter  ho;^^ 
small  the  m.esh  width  m.ight  be  taken,  since  ho  ignores  the 
facts  pointed  out  here.} 
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We  shall  introduce  tvio   different  ■schemes  employing  the 
raethod  of  finite  differences  in  a  fixed  rectangular  net  of  the 
::,t-plane.   The  first  of  those  rr.akes  use  of  the  differential 
equations  in  the  form  given  oj    O*?)?  and  we  no  lon;;cr  sripposo 
that  the  function  E  is  restricted  in  any  way.   (It  inlght  be 
noted  that  the  slopes  of  t]io    characteristics  as  given  by  (3'<3) 
are  determined  by  the  quantities  v±  c,  no  matter  iiov;  the 
function  E  is  defined,  and  hence  we  are  in  o.   position  to 
determine  appropriate  lenr^-ths  for  the  t-intervals,  in  accord 
with  the  above  discussioUj  in  the  most  general  case.   This  is 
also  a  good  reason  for  X\forl':ing  with  the  quantity  c  in  place 
of  y.)   At  the  same  tlmc^  the  calculation  is  based  on  ass-uming 
that  the  approximate  valvies  of  c  and  v  have  boon  calculated 
at  the  net  points  L,  i':,  R  (  cf .  Pig.  L!-.2)  and  that  the 
differential  equations  arc  to  be  used  to  advance  the 
A  t- 


L 


T> 


M 


/St 


-  AA 


'-i — A\ 


Pig.  k-2. 
A  Rectangular  Not 


'   « 


28. 

approximate  solution   to  the  point  P.   The  differential 
equations  to  bo  solved  are  thus 

(!:..l)         2{(c+v)c,^-:-  c^J  -;-  {(c+v)v^+  vj  +  E  =  0 

(1l.2)        -2[(-c+v)c^..  -fc^'^  +  5.(-c+v)v,,+  v^^  +  E   ^-   0 

and  the  characteristic  directions  are  determined  by  ^  =  v±  c. 

The  characteristic  vjith  slope  v+c  we  cell  the  forward. 

characteristic,  and  that  v;ith  slope  v-c  the  baclcward 

characterist5-c.   We  shall  replace  the  derivatives  in  the 

equations  by  difforcnce  o/aotients  which  approximate  the  values 

or  the  derivatives  at  the  point  11,      In  order  to  advance  the 

values  of  v  and  c  from  the  points  L,  M,  R  to  the  point  p  by 

using  (U.l)  -nd  (4.2)  it  is  natural  to  replace  the  tine 

derivatives  v,  and  c,  by  the  followino;  difference  ouoticnts: 
t      t   "  -^ 


(l;-.3)      v^ 


IlIZi^  -  ^p"  ""y. 

At  '     '       ^t  =    At 


in  both  equations.   However,  in  order  to  Insure  the  convergence 
of  the  approximations  to  the  exact  soli^.tion  when  Ax  -^   0  and 
At  — >   0  (see  [2]  for  a  proof  of  this  fact)  it  is  necessary 
to  replace  the  derivatives  v   and  c   by  difference  quotients 

XX 

vrhich   are    defined   dif f oreiitly   for    (i;.l)    than   for    ('';..'-)  ^    as 

follows: 

V     -  V  c     -  c_ 

v„  -  V,,       ■  c_  -  C-^ 
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The  reason  for    thif?  procedure  is,  c.t  bottom,  th.\t  (Ii.l}  is  an 
v^qiiation  associated  witli  tlic  forward  characteristic,  while 
([|..2)  is  associated  x/ith  the  bacI^.x^7ard  characteristic.   The 
coefficients  of  the  derivatives  and  the  function  E  arc,  of 
covxrso,  evaluated  at  the  point  M.   The  difference  equations 
replacing  (k.l)  and  (/.l.S)  are  thus  given  by 


+  4(c,.+v,.J-4^^+-4-—  r+  2(v,,,c--)  =  0 


c„-c,-   c„-c. 


Wo  observe  that  the  txvo  unlinowns,  Vp  and  c-,,  occur  linearly 
in  these  equations;  honce  they  arc  easily  found  by  solving  the 
equations.   The  result  is 

(k.8)    Vp  =  V^^  +  ^[(C^^4-V..)(|  Vj.-|  v,^-f  C-^-C^^) 

(11.9)   Cp  =  c,,+  ^  ^[{c.^^^v^^{^   v^-|  v^  +  Cj^-c^j) 
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Ill  accordance  wltli   tlao   rvTiiarks  made    above,    vio  rnust   aloo 
roquiro   that    the   r:..tio   oT   At   to  (Ax   be    taken   sraall   onon.gh 
so   that   F    lies  witlixn  tlir,    trianglo   rorracd  by  dravjlii,-/   linos 
from  L  and  R    in   the    directions    of   the    forward   and  backward 
chax'actGristics   ro::pcGtively,    i.e.    linos   with   the    slopes 
V   -*-c,    at   L  and  v   -c      at  R:      A  condition  that   is  x-foll- 

L       J..-  n       ri 

dotcrrninod  since  the  values  of  v  and  c  are  presumably  known 
at  r,  and  R . 

One  can  now  see  in  rionoral  terms  how  the  inii  ial  value 
I'^roblom  startinp-  at  t  -    '}   can.  be  solved  approximately:   One 
starts  with  a  net  alona  the  :c-axis  with  spacing  ■A'.*      Since 
c  and  V  are  knovm  at  all  of  those  points,  the  v^iluo::  of  c 
and  V  can  be  advanced  thrcufrh  use  of  (a.6)  --.nd  (i-'-.7)  to  a 
parallel  row  of  pointe  alon^  a  lino  distant  .^t  alon---';  the 
t-axis  from  the  x-axis,   'Towevor,  the  mesh  width  At  Must  bo 
chosen  small  enough  so  that  the  convorgonco  condition  is 
satisfied  at  all  net  point:;  whore  new  valuos  of  v  and  c  are 
computed. 

V/e  can  nov;  see  also  how  to  take  care  of  boundary 
conelitions,  i.e.  of  conditions  imposed  at  a  fixed  point  (say 
at  the  origin,  x  =  0)  as  given  fi.5-nctions  of  the  time.   For 
oxaraplG,  the  depth  y  (corresponding  to  the  stage  of  the  river) 
or  the  velocity  v  (which  together  x^rith  the  cross-section  area 
A  fixes  the  rate  of  discharge)  might  be  given  in  terras  of  the 
time.   Initial  conditions  dovmstream  from  this  looint  (i.e.  for 
X  >  0)  might  also  be  proscribed.   Suppose,  for  example,  that 
the  stage  of  the  river  is  pi-escribcd  at  x  =  0,  i.e.  that  y(0,t) 
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ic  known,  and  that  th.,  c-.lculation  h:-.d  already  ■•-irogrossod  30 
far  as  to  yield  values  of  v  and  c  at  not  points  alon^j  a 
certain  lino  parallel  to  tho  x-axis  and  containing  the  points 
L,  H,  R,  as  indicated  in  Fig.  ij..3.   It  is  clear  that  the 
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Fig.  U-.3 
Satisfy i n •';  E oundary  Condi^t ions 

do termination  of  the  vainer  of  v  and  c  at  point  F  can  be 
obtained  from  their  valiios  at  L,  M,  R  by  using  (ii.J)  and 
(/.;., 9)  s  as  in  the  above  discussion  of  tho  inibial  valixo  problem, 
and  similarly  for  points  ?-,  ,  Po,  etc.   Hovrcvor,  the  value  of 
V  at  Q  must  be  dctormip.od  in  a  different  manner;  for  this 
pui-posc  V7C  use  the  equation  ('''.?)  x-zith  the  subscript  Q 
replacing  P,  L  replacing  I'l,  and  M  replacing  R.   Since  v„,  c^, 
v^ ,  c„  are  supposed  knoim,  and  c^  is  also  known  since  the 
values  of  y  are  proscribed  on  the  t-axis,  it  follow?  that 
ccmaticn  (k.?)  contains  v,.  as  the  onljr  unknown;  in  fcact  it  is 
given  by  the  equation 
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Tho  reason  for  using  {I., J)    Instead  of  {li,»6)    is,  of  oourse, 
that  the  points  M  and  0  arc  associat.od  with  the  bac'rviard 
characteristic,  and  hence  (l'-.2)  should  be  used  to  approximate 
the  jc-derivatives  at  point  L  (where  the  differential  equations 
are  replaced  by  difference  eqtiations).   It  is  quite  clear  that 
the  same  general  procedure  could  be  usud  to  calculate  c„  if 
the  values  of  v  had  been  assurned  given  along  the  t-a:;is.   If, 
on  tho  other  hand,  we  had  a  boundary  condition  w^n  tjic  right 
instead  of  on  the  loft,  as  above,  wo  could  make  use  of  (.'4.6) 
for  the  forward  charactoristic  as  a  basis  for-  obtaining  the 
formula  for  advancing  uhu  solution  along  the  t-axis. 

The  above  discussioii  '.:ould  seem  to  imply  that  under 
all  c ir CUTIS tanc OS  only  one  i30iuidary  condition  could  bo  imposed 
--that  is,  that  cither  v  or  c  could  bo  prescribed  at  a  fixed 
point  on  the  river,  but  not  both-- since  prescribing  one  of 
those  quantities  leads  to  a  unique  dc terriinatlon  of  the  otlaor. 
This  is,  indeed,  truo  in  any  ordinary  ri^'or,  but  not 
necessarily  in  all  cases.   In  factj  v;c  made  a  tacit  assiomption 
in  the  above  discussion,  i.e.  that  of  the  two  characteristics 
issuing  from  any  point  of  the  t-axis  only  the  £ovyj-:jrd. 
charactoristic  goes  into  the  region  x  >  0  to  uie  rlgiit  of  the 
t-axls,  and  this  in  turn   implies  that  v+c  and  v-c,  which  fix 
the  slopes  of  the  characteristics,  are  opposite  in  sign.   The 
ph-"-sical  interpretation  of  this  is  that  the  value  of  v  (v/hich 
is  positive  hero)  inust  be  less  than  c  =  Jgj,    i.e.  th.-.t  the  flow 
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r.iust  bo  what  is  called  tranquil,  or  sub-critical/  '  Other- 
wise, as  wo  see  from  Pig.  ^-.i4-,  wc  should  expect  to  determine 
the  values  of  v  and  c  at  points  close  to  and  to  the  right  of 
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Fig,  k.li 
A  Case  of  S up c r -critical  F 1 ow 

the  t-axisj  say  at  ?:,  oj   utilizing  values  for  both  v  and  c 
along  the  segment  LQ.,  its  domain  of  dependence.   The  scheme 
outlined  above  would  therefore  have  to  bo  modified  under  such 
circiunstanccs .   Later  on,  vjo  may  wish  to  treat  flows  in 
stroam.s  so  precipitous  that  super-critical  velocities  arise 
(as  they  do  in  the  upper  reaches  of  the  New  River  In   V'est 
Virginia,  and  the-  Los  Angeles  River  even  in  its  lower  section) 
and  appropriate  finite  difference  schemes  for  such  cases  will 
then  bo  introduced.   One  sees,  however,  how  useful  the  theorjr 
based  on  the  characteristics  can  be  even  though  no  direct  use 




In  gas  dynamics  the  flow  in  an  analogous  case  would  be 
called  sub-sonic. 
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of  it  is  made  in  tho  niuncrical  calculations  (aside  frora 
decisions  on  tho  size  of  tho  t-interval). 

The  procedure  sketched  out  above,  iijhilc  it  is  pocortimendcd 
for  use  when  a  boujidary  condition  is  to  be  satisfied,  is  not 
the  best  one  to  use  for  advancing  tho  solution  to  such  points 
as  P,P  ,Pp,..,  in  Pig.  '.i..3.   For  such  "interior  points"  a 
staggered  rectangular  net,  as  indicated  in  Fig.  ii-.5:>  ap„d  a 
difference  equation  scheme  Leased  on  tho  original  differential 
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A  Staggered  Net 


equations  (3.6)  is  preferable  (cf.  [3]  for  a  discussion  of 
this  scheme).   The  equations  (3.6)  were 


(If. 11) 


2c  c..  +v,.  +  vv..  +  E  =  0    J 


cv  +  2vc  +  2c^  =  0 
JC      X     t 
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The  values  v..  and   c,.   .'t   i;ho   mid-point  M   (which   is,   ho'i/over, 
not   a  not   point)    of   the    sogriicnt   LH   arc   defined  by   the   averages: 


(1.1-.12) 
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after  which  the  derivatives  at  M  are  approximated  in  a  quite 
natuxal  way,  by  the  difference  quotients 


(^-.l3) 


V  -  v 
R    L 
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At 


Upon  substitution  of  those  quantities  into  (i.i.11).  evaluation 
of  the  coefficients  c,  v,  and  E  at  point  M,  and  subsequent 
solution  of  the  two  equations  for  Vp  and  c_,  the  result  is 


Ax 


(ii.iiij 


At,,, 


C-.  =  Cj„.  +  ^  ~-L  r_  I  C^   -  Cp^  )  V,^  +  (  V^  -  V^  )  C.,-1 


As  we  see  on  comparison  uith  (l|-.8)  and  {.l^-.g/j  these  equations 
are  simpler  than  the  earlier  ones.   The  criterion  for 
convergence  remains  the  sa;-i.e  as  before j  i.e.  that  ?  should  lie 
within  a  triangle  formed  by  the  segment  LR  and  the  two 


charactor'i sties  issuing  frciii  its  ends. 
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